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Abstract 

The paper studies some inverse boundary value problem for simplest parabolic 
equations such that the homogenuous Cauchy condition is ill posed at initial time. 
Some regularity of the solution is established for a wide class of boundary value 
inputs. 
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Parabolic equations such as heat equations have fundamental significance for natural 
sciences, and various boundary value problems for them were widely studied including 
well-posed problems as well as the so-called inverse and ill-posed problems that often are 
very significant for applications (see, e.g., Beck (1985)). The present paper investigates 
a inverse boundary value problem on semi-plane for homogenuous parabolic equations 
with homogenuous Cauchy condition at initial time and with Dirichlet condition on the 
boundary of semi-plane. The parabolic equation is the equation of a backward type 
that is usually solvable with Cauchy condition at terminal time. However, we consider 
this equation with Cauchy condition at initial time so that the problem is an inverse 
problem. One may think that these problems are always ill-posed in the sense that 
there is no regularity of solutions such as prior estimates for the solution via some 
norms of free terms (see, e.g., Beck (1985), Tikhonov and Arsenin (1977)). We found 
a class of very generic parabolic equation with constant coefficients and with certain 
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sign of the drift coefficient (i.e., the coefficient for first derivative) such that the inverse 
problem has some regularity. More precisely, we found a wide enough class of inputs in 
the boundary value condition on the boundary of the semi-plane that ensures regularity 
in a form of prior energy type estimates. This class of inputs is everywhere dense in 
the class of L2-integrable functions; it includes differentiable functions. 



1 The problem setting 



Let us consider the following boundary value problem on semi-plane: 

du d 2 u du 

a— (a?, t) + — — k (x, t) + b—(x, t) + cu(x, t) = 0, x > 0, t > 0, 
at ox z ox 

u(x, 0) = 0, x > 0, 

du 

k u(0,t)+k 1 — {0,t)=g(t), t>0. 
ox 

Here x > 0, t > 0, and a, b, c,ko,k\ G R are constants. 
We assume that 



(1) 



a > 0, b > 0, kl + h\ > 0, fc &i < 0. 



(2) 



The assumption that a > and the presence of the initial condition at t = makes 
problem (pQ) an inverse problem (see, e.g., Beck (1985), Tikhonov and Arsenin (1977)). 
We assume that u(x, t) = and g(t) = for t < 0. 

Let r denotes the set of all functions g : R — > R such that g(i) = for t < and 
with finite norm 

<9g 



ISllwlfR) — H5lli 2 (R) + 



dt 



L 2 (R) 



Remark 1 Functions g £ T are continuous on R and vanising on t < since 
dg(t)/dt £ 1,2 (R). -For instance, g(t) = e~ l smt £ T, but g(t) = e~ l cost ^ T. 



Let fl = Rx R + . Let W be the space of the functions v = v(x, t) : R X R -1 



R 



such that v(x,t) =0 for t < and with finite norm 



v w 



Il 2 (d) + 



<9it 
9a; 



+ 



d 2 v 
dx 2 



L 2 (D) 



+ 



dv 

di 



L 2 (D) 



The class W is such that all the equations presented in problem ([T]) are well defined 
for any u G W. Let us show this. If v G W, then, for any i* > 0, we have that 
^Ir+xIo,*,] ^ C([0, t*], L2(R + )) as a function of t G [0, t*]. Hence the initial condition 
at time t = is well defined as an equality in L2(R + ). Further, for any x* > 0, we have 
that f|[o,a;«]xR+ G C([0> x*], L 2 (R + )) and || G C([0, x*], L 2 (R + )) as functions 



of x G [0,x*]. Hence the functions go(i) = u(0, £), <?i(t) = ^(x, £)| x= o are well defined 



as elements of L2(R + ), and the boundary value condition at x = is well defined as 
an equality in L2(R + ). 



2 The main result 

Theorem 1 Let condition flj) be satisfied, and let 

H = b 2 /4 < c. (3) 

Then there exists a unique solution u(x,t) in the class W of problem {7]j in the domain 
D for any g G T. Moreover, there exists a constant C = C(a,b,c,ko,ki) such that 

\\u\\w < C\\g\\ w i (R) . (4) 

Remark 2 It will be seen from the proof that it is crucial that u(x, 0) = in the 
Cauchy condition and that the parabolic equation is homogenuous. We cannot extend 
the result for non-zero initial conditions or non-zero free term in the parabolic equation. 

The following theorem shows that assumption ([3]) is not really restrictive if we are 
not interested in the properties of the solutions for T — > +oo, as can happen if we deal 
with solutions a finite time interval. 

Theorem 2 Let condition (0j holds, but (0j does not hold. Let M be such that 
g(t)e~ Mt G r and 6 2 /4 < c + M. Then problem ^) has a unique solution u such 
that um G M., where uu{x,t) = e~ Mt u(x,t). 

Proof of Theorem^ Let R+ = [0, +oo), C+ = {z G C : Re z > 0}. For v G L 2 (R), 
we denote by Tv and Cv the Fourier and the Laplace transforms respectively 

V(iLo) = {Tv){iuj) = I e~ iujt v(t)dt, to G R, (5) 
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V(p) = (Cv)(p) 



i r°° 

= -4= / e- pt v(t)dt, p E C+. 



(6) 



Let H r be the Hardy space of holomorphic on C + functions h(p) with finite norm 
\ H r = sup fc>0 \\h(k + iu})\\ L r(^, r E [1,= oo] (see, e.g., Duren (1970)). 
Let u E W be a solution of Q. Set g (t) = u(0,t), gi(t) = ^(x,t)\ x=Q . As was 
discussed above, the functions g\. are well defined as elements of L2(R + ). 

Let G = Cg. We have that G E -ff 2 . Let Gfc(p) = and J7 = £u be defined for 
p E C + . They are well defined since u E W; in addition, E -ff 2 . 

For functions V : R + x C + — ► C, where C + = {z : Rez > 0}, we introduce 
semi-norms 



\V\ 



\V\ 



L22 



w* 



R 



dx J \V(x,iu>\ du 



1/2 



\V\ 



dV 



dx 



+ 



d 2 v 



and norms 



IVI 



L 22 



R+ 



\\v{x,-)r H 2dx\ 



dx 2 

s 1/2 



£22 



+ 



SF 
' at" 



I^IIh = II^IIlh 



Instead of (JT]), consider the problem 



dV 



dx 



1 11 

L 22 



+ 



d 2 v 



dx 2 



^22 



ay 



at 



^22 



d 2 U dU 
apU(x, p) + — —~ (x, p) + 0-^— (a;, p) + cU(x, t) = 0, x > 0, 
ox z ox 

BIT 

k U(0,p) + —(0,p) = G(p), p E C + 



subject to the following condition 
Sf/ 5 2 ?7 



|C/||-H < +00. 



(7) 



(8) 



Let Afc = Afc(p) be the roots of the equation A 2 + 6 A + (c + ap) = defined for 
p E C + as Ai = —6/2 — y/Jl — ap and A2 = —6/2 + ^//x — ap, where /i = 6 2 /4 — c < 0. 
We mean the branch of the square root such that Arg yfjx — ap E [— 7r/2, +vt/2] and 
Re ^//x — ap > 0. Under these assumptions, the function -^//x — ap is holomorphic and 
does not have zeros in C + . We have that 



ReAi(p)<--, pEC+, 

35 > 0, > : ReA 2 (iu;) > 5 if u E R, > 



0) 
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In addition, we have that the functions Xk{p) are holomorphic in C + , and 



(Aib) - Hp))' 1 g h™, Aft(p)(Ai(p) - \ 2 (p))- 1 G fe = 1, 2, 

(fco + fciAiCp))- 1 g#°°, Aifo + fciAiCp)) -1 eff 00 . (10) 



The last two statements here follow from (121. Let 



iV 



Ai — Ao 



jfc=l,2 



Ai — A2 



+ 



1 



k + /ciAi 



+ 



Ai 



fc + feiAi 



#00 



It can be seen also that the functions e xXk ^ are holomorphic in C + for any x > 0. 
For any x > 0, the unique solution of ([7]) is 

1 



U(x,p) 



(G^p) - X 2 G (p))e x ^ + (Gi(p) - A!Go(p))e A2:E ) . 



A 2 x 



Ai — A2 

This can be derived, for instance, using Laplace transform method applied to linear 
ordinary differential equation (J7]), and having in mind that 

1 1 1/11 



X 2 + b\ + c-ap (A-Ai)(A-A 2 ) A1-A2V.A-A1 A - A 2 / ' 
A A 1 ( Ai A 2 



A 2 + 6A + c-ap (A-Ai)(A-A 2 ) Ai-A 2 \A-Ai A - A 2 
Let represent U as U(x,p) = U\(x,p) + U 2 (x,p), where 



U 1 (x,p)=e x ' x J 1 (p), ,h{p) 
U 2 (x,p)=e x * x J 2 (p), J 2 (p) 



1 



Ai — A2 
1 



(Gi(p)-A 2 G (p)), 
(Gi(p)-AiG„(p)). 



Ai — A2 

By ([9]), the fact that u G W implies that ([8]) holds. Let us show that 

Il^l||i22 < +°°- 

By ®, |e xAl (P))| < e"^/ 2 < 1, p G C+ It follows that 

||e xAl(p) Ji(p)|| H 2 < sup le^^HlJiUfla < e- fex/2 ||Ji|| H 2 <7V e -^/ 2 ^ ||G fe ||^2. 

P eC+ fc=0,l 



(11) 



Then ® and ^ imply that 



|^2||l 2 2 < +°°- 



(12) 



Further, 



R+ 



dx f \U 2 {x,iuj)\ 2 duj= f dx [ \ e x ^ iuj > ' J 2 {iu)\ 2 'du 



dx e RcX2{luj)x \J 2 {iuj)\ 2 duj> / dxe 5x / \J 2 {iuj)\ 2 du . 

Note that the J 2 (iio) is vanishing on {uj : \uj\ > uj*}. Since J 2 G -ff 2 , it follows that 

«/ 2 = (G 1 (p)-A 1 G (p))e A2:E = 0, 



i.e. 



Gx(p) = XiG (p). 
Remind that koGo(p) + kiG±(p) = G{p). It follows that 

k G (p) + ki\xG (p) = G(p), G (p) = (ko + hX^p^Gip), 
hip) = G ( P ), U(x,p) = Ufap) = e x ' x Go(p). 



Let us estimate 

By @, \e xXl( -P»\ < e~ bx / 2 < 1. It follows that 



\\p m e xX ^G (p)\\ H2 < e~ bx / 2 \\p m G (p)\\ H2 

<e- bx / 2 \\(k + k 1 X 1 (p))- 1 \\ H o O \\p m G\\ H 2 <Ne- bx / 2 \\p m G\\ H i, m = 0, 1. 



It follows from the above estimate that 



\\p m U\\ L n ■<NC 1 \\p m G\\ m , m = 0,l. 



Further, we have that 



A 



r)TI 

^L( x , p )=G (p)X 1 e x ^ x = , 

OX Kq + k\X\ 



1 -—G{p)X 1 e XlX . 



We obtain again that 



dU 



dx 



L?2 4 



dx < NC 2 
H 2 -'R 4 



,\i{p)x 



G( P ) 



H 2 



dx 



< NC 2 



R 4 



- bx / 2 \\G{p)\\ 2 H2 dx<Cz\\G\\ 



H 2 • 



(13) 



(14) 



(15) 
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By ((7D, d 2 U/dx 2 can be expressed as a linear combination of U, pU, and dU/dx. By 

(HHD-dlSD, 



7r+ <9x 2 
It follows that 



dx < C4 



H2 




017 

dx 



(x,p) 



dx+ ^2 

m=0,l ^ R4 



|p m tf(s,p)||^dx 



R+ 



<9x 2 



(x,p) 



ff2 



^<C 5 (||G|| 2 H2 + |bG(p)||^). 



(16) 



Here Ck are constants that depend on a,b,c,ko,ki. By (fT3|) - (fT6|) . estimate flSJ) holds. 

Let u(x,-) = ^ r_1 C/(x, iw)| a , e R. By fp~3|) . it follows that the corresponding inverse 
Fourier transforms ti(x, •) = J : ^ 1 U{x, iw)| wg R, §r(av) = F~ l {pU(x, iu)\ u gR.) are well 



defined and are vanishing for f < 0. In addition, we have that U(x, iuj) = U(x, —ioj) 
(for instance, Gk(iuj) = Gk(-iuj), e xX *( iu} ) = e xX k(-^) ^ e ^ c ^ ^ f u ows that the inverse 
of Fourier transform u(x, •) = J-~ 1 U(x, •) is real. By ([8]), estimate @ holds. Therefore, 
u is the solution of (H|) in W. The uniqueness is ensured by the linearity of the problem, 
by estimate (j4|), and by the fact that Cu(x, •), C(d k u(x, -)/dx k ), and C(du(x, ■ /dt) are 
well defined on C + for any function u from W, i.e, ([7|) must be satisfied together with 
([8]) . This completes the proof of Theorem [TJ □ 

Proof of Theorem [|J Rewrite the parabolic equation as the one with c replaced by 
c + M and g(t) replaced by g(t)e~ Mt . By Theorem [IJ solution um £ .M of the new 
equation exists. Clearly, u(x, t) = e Mt UM(x, t) is the solution of of the original problem. 
□ 



3 Some application 

Total absorbing on the boundary 

Let T > be given. Let us consider the following well-posed boundary value problem 
on semi-plane: 



dv 



d 2 v 



. . dv 

— (x, t) + a— -5 (x, t) + b— (x, t) + cv(x, t)=0, x > 0, t £ [0, T], 
at ox z ox 

v(x, T) = v*(x), 

dv 

k v(0,t)+k l — (0,t) = g(t). 



(17) 



Here a, b, c,ko,ki G R are constants such that holds. 

Theorem 3 For any g G T , i/iere exists G L2(R + ) smc/i i/iai i>(:r,0) = ; where 



v G W is t/ie solution of well-posed problem JT7[), 



Proof. It suffices to take the solution u G W of problem ([!]) and take = u(x,T), 

v = u. □ 

Note that one can rewrite problem (I17D as a well posed problem for forward parabolic 
equation with initial time at time t = via time change t —* T — t. In that case, the 
phenomena described in Theorem [3] looks more impressive. 

Restoring past distributions of diffusion processes 

Consider the following stochastic process 

y x (t) =x + bt + aw(t). (18) 

Here x > 0, w(t) is a scalar Wiener process, b > and cr > are constants. 

Let a be a random number such that a > and it has the probability density 
function p G L2(R + ) which is supposed to be unknown. We assume also that a is 
independent from w(t) — w(t\) for all t > t\ > 0. 

Let y(t) = y a (t) be the solution of Ito equation (|18f) with the initial condition 
y(0) = a. Set r a = min{t > : y a (t) = 0}. 

Let I denotes the indicator function of an event. 

Let p(x,t) be the probability density function of the process y a {t) if this process is 
killed at and inside (0, +oo) with the rate of killing c (case of c > is not excluded 
for the sake of generality) . More precisely, p is such that 



/ p(x,t)dx = Ee c \ y aM eB )Ii T a> t \ 
Jb 



for any domain B C R + . It is known that evolution of p is described by the parabolic 
equation being ajoint to ([T]) with the boundary value conditions p(x, 0) = p(x), p(0, t) 
0. 

For g G T, let u g = u g (x, t) be the solution of the problem ([T]), where ko = 1, k± = 0. 
Let T > be given, and let *& g (x) = u g (x,T). 
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Theorem 4 For all functions g G T and all c E R, 

Ee CTa 5 (^)V<T} = -Ee cT ^(y(T))I {Ta > T} = - / p(x, T)* g (x)dx. (19) 

Theorem fallows to solve effectively the following inverse problem: find the distribution 
of T a I{ T a <T j for unknown distribution of a using the "future" values of p(x,T) only. 
More precisely, one can find the values of the expectations at the left hand side of (jl9h 
for all g G T using only p(x, T) via the following algorithm: 

(a) Find the solution of ([!]); 

(b) Find V g = Vg(;T); 

(c) Using known p(x,T), calculate the integral at the right hand side of (fT9l) . 

This approach does not require regularization of a ill-posed problem such as in Beck 
(1985) or Tikhonov and Arsenin (1977). 
Proof of Theorem [|| Clearly, 

Ee^ a ^u g (y(^ A T), r a A T) = Ee CT " g(r a )I {ra<T} + Ee cT * 9 (y(T)) I {t „> t} . 
By Ito formula, 

Be c( - raAT ^u g {y(T a A T), r a A T) = Bu g (a, 0) = 0. 

In addition, 

Ee cT ^ 9 (y(T))I {r0 > T} = f p(x,T)V g {x)dx. 
Then the result follows. □ 
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